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ABSTRACT 

We performed simple spot-model calculations for quasi-periodic brightness 
variations of solar-type stars showing superflares, by using Kepler photometric 
data. Most of superflare stars show quasi-periodic brightness modulations with 
the typical period of one to a few tens of days. Our results indicate that these 
brightness variations of superflare stars can be explained by the rotation of the 
star with fairly large starspots. Using the result of the period analysis, we in- 
vestigated the relation between the energy and frequency of superflares and the 
rotation period. Stars with relatively slower rotation rates can still produce flares 
that are as energetic as those of more rapidly rotating stars, although the average 
flare frequency is lower for more slowly rotating stars. We found that the energy 
of superflares are related to the total coverage of starspots. The correlation be- 
tween the spot coverage and the flare energy in superflares is similar to that in 
solar flares. These results suggest that the energy of superflares can be explained 
by the magnetic energy stored around starspots. 

Subject headings: stars: activity — stars: flare — stars: rotation — stars: solar-type 
- stars: spots 
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Introduction 



Solar flares are the most energetic explosi ons on the surface of the Sun, and are thought to 
occur by release of magnetic energy (e.g., 



Shibata &: Magara 



2011). F 



ares are also known 



1989 



Gershberg 



to oc cur on various types of stars including solar-type stars ( ISchaeferl 

20051 ). Among them, young stars or close binary stars sometimes produce "superflares", 

flares whose total energy is 1 ~ 10 6 times more e nergetic (~ iq 33 ~ 38 erg) than the largest 



flares on the Sun (~ 10 erg) (jSchaefer et al. 



20001 ). Such stars generally rotate fast (v TOt ~a 



few 10 km s 1 ), and the magnetic fields of a 



large regions on the stellar surface (JGershberg 



'ew kG are considered to be distributed in 



2005 



Shibata fc Yokoyama 



1999 . 



20021 ). In 



contrast, the Sun slowly rotates (v rot ~2 km s _1 ), and the magnetic fields are weak. Here 
we define "Sun-like" stars as solar-type stars slowly rotating and whose surface tempera- 
ture is 5,600K ~ 6,000K. It has been thought that superflares cannot occur on Sun-like stars. 



Schaefer et al.l ( 120001 ). however, found 9 candidates of superflares on slowly rotating 
stars like the Sun. This is extremely important in many fields including magnetic activity 
research in solar /s tellar physics as well as the planetary habitability in astrobiology (e.g., 



Segura et al. 



20101 ). The frequency, detailed properties, and mechanism of the superflares 



are, however, still not clear because of lack of observations. Therefore, it is necessary to 
investigate in detail how often superflares occur on solar-type stars, properties of superflares, 
and stellar conditions which lead to superflares. 



We have already analyzed the data by the Kepler spacecraft ( IKoch et al. 



20101), and 



lave surface tempera ture of 



20121 ). The 



discovered 365 superflare events on 148 solar-type stars that 

5,100K < T e ff < 6,000K, and surface gravity of log g > 4.0 ( iMaehara et al. 

Kepler spacecraft is very useful to detect faint brightness increases in the stellar brightness 

due to stellar flares because Kepler realized high photometric precision and continuous 
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time-series data of a lot of stars over a long period (jWalkowicz et al.ll201ll ; iBalonal 120121 ) . 
We found that superflares whose energy is 10 2 — 10 3 times larger than that of the most 
energetic flare on the Sun, can occur on solar-type stars once in a few thousands of years. 
Many of solar-type stars having superflares show quasi-periodic brightness variations with 
the typical period from one day to a few tens of days. Such variations can be explained by 



the rotation of the star (e.g.. lBasri et al. 



2011 



Debosscher et al.l 



2011 



Harrison et al. 



2012|). 



In this paper, we investigate the brightness variations of superflare-generating stars 
by considering the signature of stellar rotation and starspots. First we show that the 
brightness variation of some typical superflare-generating stars are explained by the 
rotation of a star with starspots, by simple model analyses. Our main purpose of the 
model analysis is to demonstrate that the brightness variations include the information 
of the overall areal coverage by starspots and the rotation periods of superflare stars. 
Second, using the result of the period analysis, we consider relations between the rotation 
period and properties of superflares such as flare energy and flare frequency. Third, 
we discuss relations between the starspot coverage and the energy of superflares by 
assuming that the brightness variations are due to the rotation. For the Sun, it has 
been known that there is a positive correlation between the sunspot coverage and the 



energy of the largest flare observed, and that the energy of the largest : 



Sammis et al. 



are d epends 



20001 ). We 



on the magnetic energy stored around the sunspots (e.g. 

then investigate whether this correlation can be applied to superflare-generating stars, 

and discuss the relations between the magnetic energy of starspots and the superflare energy. 



The details of Kepler data we use in this paper are explained in Section [2j We 
show in Section [3] the results of the spot modeling for some typical superflare-generating 
stars. We consider the relation between the rotation period and properties of superflares, 
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on the basis of the results of the period analysis in Section HJ We discuss the relation 
between the starspot coverage and the energy of superflares by assuming that the brightness 
variations of superflare-generating stars are due to the rotation in Section 

2. Observational Data and Period Analysis 

We searched for flares on solar-type stars using the Kepler data which was taken during 
the period from April 2009 to September 2010 (quarter 0~6). These data were retrieved 
from the Multimission Archive at Space Telescope Science Institute (MAST)Q|. We used 
the effective tempe rature (T eS ) and the surface gravity (logg) in the Kepler Input Catalog 



( JBrown et al.l 120111 ) to select solar-type stars. The selection criteria are as follows; 
5, 100 < T e ff < 6, 000, and logg > 4.0. The total number of solar-type stars is about 90,000. 
The length of the observation period during each quarter and the number of solar-type 
stars are summarized in Table [TJ 

We analyzed the lon g-cadence (time res olution of about 30 min) flux detrended by 



PDC-MAP pipeline ( IStumpe et al 



detection method are described in 



2012 }) for det e ction of su perflares. Details on th e 



Maehara et al. 



(120121 ) and 



Shibayama et al. 



(120131 ) . The 



version of the data we use in this paper are shown in Table [Q 

Figure Q] shows the typical light curves of solar-type stars having superflares (KIC6034120 
& KIC6691930, Quarter 2; KIC10528093, Quarter 0-2). KIC6034120 shows a lightcurve 
of a simple sine-like profile which can be well reproduced with one large starspot. The 
lightcurve of KIC6691930 requires two starspots to explain the shoulder of the rising 



^ttp:/ /archive. stsci.edu/kepler/ 
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part. KIC10528093 also needs two starspots, and its long-term amplitude variation cannot 
be explained by the solid rotation model. These parts of the data are selected to make 
these featur es clearly demons t rated . Light curves of all solar-type stars having superflares 



reported in 



Shibayama et al. 



(120131 ) are shown in Online-only Figures. The period of the 



brightness variation was estimated by the discrete Fourier transform (DFT) method. We 
carried out the following 5 procedures: (1) We measured the standard deviations (er) of the 
data in each quarter, and deleted all the data of a quarter whose a value is over three times 
larger than the mean value. (2) The linear trend was removed in the brightness during each 
quarter. (3) The gaps in the mean brightness between each quarter were adjusted. (4) We 
took the power spectra of the period from 0.1 day to 45 day, since the reduction pipeline 
seems to generate spurious peaks at timescales longer than 45 days, a half of typical lengths 



of each quarter 



(5) We chose the peak in the power spectrum whose amplitude has the 



highest ratio to the red noise spectrum (e.g.. 



Press 



1978 



Vaughan 



20051 ) as the period of the 



brightness variation of the star, since in some cases simply taking the highest peak can lead 
to choose spurious peaks at long timescales. Figure [2] is the results of the period analysis 
for the brightness variation shown in Figure [TJ Such stars show quasi-periodic brightness 
variations with the typical period from one day to a few tens of days. 



2 For the detrend method by the Kepler pipeline, see 



http://keplergo.arc.nasa.gov/PyKEprimerCBVs.shtml 
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3. Model Calculation of Brightness Variation of Rotating Stars with Large 

Starspots 

We performed simple calculations to show the example that brightness variations of the 
superflare-generating stars can be explained with the rotation of these stars having large 
starspots. There are many former researches which calculate the br ightness variat i on by 



assuming the exi stence of dark spots on the surface of the star (e.g.. 



1987 



Eker 



Budding 



1977 



Dorr en 



19941 ). Spot modeling has much uncer tainty and i t is almost impossible to 



get unique solutions fro m the spot modeling (e.g., 



Walkowicz et al. 



Eker] (119961 ): 



Kovari fc Bartusl (119971 ): 



(120131 )). There are many parameters such as inclination angle and spot 
latitude which affect the shape of the lightcurve. In addition, the Kepler spacecraft provides 
only single color photometry data, and this also generates degeneracy between the spot size 
and the spot contrast. Our main purpose is to demonstrate that the brightness variations 
include the information of the overall areal coverage by starspots and the rotation periods 
of superflare stars. 



Our method of the spot modeling is as follows. We regard the total luminosity of 
the photosphere, and the starspot area as a good representation of the flux (F) in our 
modeling. The shape of the model star is assumed to be a sphere, and some circular spots 
are placed on the photosphere. It should be noted that we do not exclude the possibility 
that starspots of superflare stars are collections of smaller spots. The surface of the star is 
divided into 180 pieces (J) in the latitudinal angle and into 360 pieces (k) in the longitudinal 
angle. The contributions to the total brightness are summed up from each piece {fj t k) as 
follows; 

F<*T,fi*, (1) 



and 



f, 



f s ^, ot (in the spot area) 



3,k 



pphot 



ffk° ( ou ^ °f the spot area) 



(2) 



The temperature of the photosphere (T p hot) and that of the spot area (T spot ) are here 
considered t o be the same oy er the solar photosphere (6,000K) and the sunspot area 



(4,000K, cf. 



Berdyugina 






(3) 



20051 ) , respectively. The luminosity of each piece which is in the 
starspot area is estimated by 

_,„,, j crsBr s 4 pot S'pi X n j>k ■ aj, k (n j>h ■ a j)k > 0) 
[ (n i)fe • a i)fe < 0) , 

where S P [ X and rij^ is a size and a normal vector of each pixel, respectively, and a^ k is a 
unit vector which is antiparallel to the line of sight. <tsb is the Stefan-Boltzmann constant. 
The luminosity of each piece which is out of the starspot area is estimated by 

cr SBr phot S , pix n jtk ■ a j)k (n j>k ■ a,j )k > 0) 
(n.j, k ■ a jjk < 0) . 



^■phot 
■>j,k 



(4) 



We also take into the effect of limb darkening. Applying the equation (1) in ISingi (120 lOl ). 
this effect is calculated by 

1 — u(l — cos/i) , (5) 



UM 



fNLD 
Jj,k 

where / ? N ^ D is the luminosity of each piece without the limb-darkening effect, fj y k{^) is that 
with the limb-darkening effect, jjl — cos 6 and 6 is the viewing angle. We here assume 



it 



0.6 on the basis of Figure 3 in 



Claret 



(120041 ) . Because of the rotation, the visibility 



of the starspots changes, and the brightness then varies. We calculate the brightness by 
changing the angle of rotation by 1 degree. Model light curves are drawn with arbitrarily 
changing the spot position, spot radius, and inclination angle i, which is the angle between 
the line of sight and the rotation axis of the star. We selected cases which show rough 
agreement between the model and the observation by eye. 
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Figures [31 0], and [5] show results of the spot model calculation for the brightness 
variations of typical solar-type stars having superflares in Figure [T] (KIC6034120, 
KIC6691930, and KIC10528093). The stellar parameters of these stars are shown in Table 
|2j The period of the brightness variations of these stars is about 5~15 days, and the 
amplitude is from 0.1% to a few %. 

As shown in Figure [31 the simple sine-like light curve of KIC6034120 can be reproduced 
by the model with one spot. Figure [3] (a) shows the model light curve for KIC6034120. 
The best set of model parameters are listed in Table |3J Snapshots of the modeled star 
are displayed in Figure[3] (b) and (c), which shows the star has a large spot compared to 
the Sun. Figure[3] (d) shows the comparison of observed light curve and the model light curve. 

The shape of the light curve of KIC6691930 has one peak and one shoulder as 
shown in Figured] (b). The features can be reproduced by the two spot model as shown in 
Figure HJ The longitude difference between two spots can be estimated by measuring the 
rotational phase of the shoulder in the light curve. Figure H] (a) shows the model light curve 
for KIC6691930. The model parameters are listed in Table |3j Figure H] (b), (c), and (d) are 
drawn for KIC6691930 in the same manner of Figure [3] (b), (c), and (d) for KIC6034120. 

The shape of the light curve of KIC10528093 also has a shoulder as shown in Fig- 
ure [H (c). However, the long-term amplitude variation cannot be reproduced by a simple 
model with two spots in solid body rotation since the phase of the shoulder changes as time 
progresses. 



Thes e variations can be explained by the differential rotation (e.g., 



Frasca et al. 



2011 



Frohlich et al. 



20121 ) . Here, we consider the two-spot model with the differential 
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rotation. It should be noted that what this model shows is only that the differential rotation 
in this model is of approximately the correct amount. This is because a qualitatively 
similar fit would be produced by any models, with almost any spot parameters, which 
get the amount of the differential rotation correct and the rough longitudinal separation 
correct. Parameters of the differential rotation (the depend ence of the rotation spee d on 



the latitude) we use here is based on the value for the Sun (JSnodgrass fc Ulrichlll990l ). The 



relation between the angular velocity of the rotation (u) and the latitude (cp) is represented 
by the equation (jSJ). 

uo = A + B sin 2 (f + C sin 4 (p . (6) 

For the Sun, the values of A, B and C are as f ollows; A = 14.71 deg d ay -1 , 



W- 



1990). We here 



B = —2.39 deg day -1 , and C = —1.78 deg day -1 (JSnodgrass fc Ulrichl 
ignore the "C sin 4 <p" term since this term is negligibly small in our model, and then we 
convert values A and B to ones which corresponds to the rotation period of KIC10528093. 
For KIC10528093, the value of A and B are as follows; A = 31.1 deg day -1 , and 
B = —5.05 deg day - . Figure [5] (a) shows the model light curve for KIC10528093. The 
model parameters are listed in Table [31 Snapshots of the modeled star are shown in Figure 
[5] (b), (c), (d), and (e). We compare the observed light curve with modeled light curve 
in Figure [5] (f). Figure [6] shows power spectra of the observed light curve and the model 
light curve in Figure [5] (f). We can see that the feature of the power spectrum, that is, 
the dominant and overtone peaks, and the slope to the higher frequencies, is also roughly 
reproduced with our spot model. 

Consequently, most of the brightness variations with the period of a few days to a 
few tens of days and the amplitude of 0.1 ~ 10% are expected to be able to be explained 
by assuming that the star has fairly large starspots. 
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4. Ensemble Properties of Superflares: Dependence of Flare Energy and Flare 

Frequency on Rotation Period 

Applying the results shown in Section [31 we here assume many of the brightness variations 
with the period of solar-type stars having superflares can be explained by the rotation of 
the star with large starspots. Only three particular lightcurves were shown to be roughly 
reproduced by our uncomplicated model calculations in the above section. This does not 
prove that such spot models can be applied to the lightcurves of all the superflare stars. 



Figure [7J shows the relations between the rotation period and the 
flares. Although the similar relations have been alr eady discussed in 



Shibayama et al. 



eatures of supe r- 



Maehara et al 



(2012) 



( 120131 ) from Kepler 



we here add data of superflares we newly found in 
quarter 3-6 data and improve the statistical precision by increasing the numbers of data. 
The rotation periods are determined by using the w ay described near the e nd of section 
El although this way includes some uncertainty (e.g.. IWalkowicz et al.ll2013l ). Figure [7J (a) 
indicates that the most energetic flare observed in a given rotation period bin does not have 
a correlation with the period of stellar rotation. If the superflare energy can be explained 
by the magnetic energy stored near the starspots (we discuss this in Section [5]), this result 
suggests that the maximum magnetic energy stored near the spot does not have a strong 
dependence on the rotation period. Figure [7J (b) shows that the average flare frequency 
in a given period bin tends to decrease as the period increases to periods longer than a 
few days. The value of the frequency is the average of all superflare stars in the same 
period bin. Some of superflare stars (e.g., examples in Figured]) are then able to have 
higher frequencies of superflares. The frequency of superflares on rapidly rotating stars is 
higher than slowly rotating stars. It is known that the rotation period correlates with the 
chro mospheric activity a nd the more rapidly ro tating stars have higher magnetic activity 



(e.g., 



Noyes et al 



1984 



Pallavicini et al. 



19811 ). These imply that rapidly rotating stars 
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with high magnetic activity can generate more frequent superflares. 



5. Discussion on starspot coverage and Energy of superflares 

As described in Section [H it has been said that the brightness variation of many solar-type 
stars observed by Kepler is due to the rotation of the star with starspots. This effect of 
rotation is also seen in the Sun. Figure [8] represents that there are also the same relations in 



the Sun; changes of the visible area of the sunspot causes the brig 
faculae also can affect the brightness variations of the Sun (e.g., 



itness variation. However, 



Lanza et al. 



20031). The 



brightness vari ations of superflare-gene rating stars are also probably somewhat affected by 



faculae, though ILockwood et al 



(120071 ) indicates that the photometric behaviors of (young) 



active stars are dominated by starspots. In th is paper, we do n o t inc 



for analyzing the Kepler data, as described in 



Walkowicz et al. 



(2013). 



ude the faculae effects 



Maehara et al. 



( 120121 ) showed that superflare-generating stars have the brightness 
variation with the period of a few days to several tens of days, and that many of such 
variations are likely due to the rotation, though the possibility that the star is a binary 
is not completely excl uded. It is well known that more than half of general solar-type 
stars are binary (e.g. 



Duquennoy fc Mayor 



19911 ). However photometric observations by 



Kepler mission cannot detect spectroscopic binaries. We cannot exclude the possibility 
that superflares are taking place on the spectroscopic binary companions. Because of this, 
the period data shown in Figure 7 are not necessarily due to rotation, and in particular 
the target with a 0.1-day periodicity is probabl y a binary. Spectro scopic observation is 



necessary, which will be our future project (e.g., 



Notsu et al. 



2013a). 



13 



Here we assume the brightness varitations are due to the rotation. Our model anal- 
yses for some typical brightness variations suggest us that, if we take into account the effects 
of some parameters, such as, the inclination angle, and the spot latitude, the brightness 
variations whose amplitude is about 0.1~10% can be well explained by the rotation of the 
star with fairly large starspots. 

In the following, we discuss the relation between the starspot coverage and the energy of 



super flares. Flares are the release of the stellar magnetic energy (e.g., 



Shibata &: Magara 



20111 ). It is known that there are posi tive correlations bet ween the sunspot coverage and 



the energy of the largest solar flares (jSammis et al. 



20001 ) . It is, therefore, important 



to consider whether the same correlations can be found and the observed maximum of 
superflare energy can be explained by the magnetic energy of large starspots which many 
of the superflare stars are expected to have. Although we have already simply discussed 
the relation between the amplitude of superflares and that o f the brightness variations in 



Supplementary Information Section of 



Maehara et al 



(J2012l ). we here advance this former 



analysis and discuss whether the magnetic energy of starspots can well explain the energy 
of superflares by applying the results in this paper that starspot coverages can be roughly 
estimated from the brightness variation. 



First, we discuss expected relations if the energy sources of superflares are the mag- 
netic energy stored around the starspots (-Efiare < -Emag)- The total energy released by the 
flare must be smaller than (or equal to) the magnetic energy stored around the starspots. 
The order of the stored magnetic energy (E m3jg ) can be roughly estimated by 

BH 3 



E « 



8^ 



(7) 
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where B and I correspond to the magnetic field strength and the size of the starspot region. 
It should be noted that we cannot exclude the possibility that starspots of superflare stars 
are collections of smaller spots. 

We roughly assume that there is a linear correlation between the brightness varia- 
tion and the spot coverage when we estimate the magnetic energy stored around the 
starspots from the amplitude of the brightness variation. The total amplitude of the 
brightness variation due to the rotation normalized by the average brightness (AF rot ) can 
be expressed by 



AF ri 



1 - 



T \ < 

- L spot \ 



.4 



spot 



.4 



phot 



- J- phot - 

where A spot is the total area of starspots, and A phot is the total visible area of the stellar 
surface. Since the Kepler photometer covers a wide spectral rage (4,000 ~ 8,500A), observed 
brightness changes would be nearly approximated by the bolometric brightness change 
obtained by the equation (jHJ). 

If we assume that / 3 « A spot , the equation (J7J) can be transformed to 



#\ 3 /2 _B 2 s A phot AF Tot \3/2, 

nla S ~ S,r s P ot ~ S^ V 1 _ (T . IT . .\4 



" tt Stt^ 87rVl-(T spot /T phot r J '~ Lll »" ' ' 



where -Eflare is the total energy released by the flare. The amplitude of the flare normalized 
by the average brightness (AFfl are ) can be estimated by 

AF flarc « ^ , (10) 

where L star is the luminosity of the star, and r is the e-folding time of the flare. Then the 
relation between the flare amplitude and the brightness variation amplitude can be written 

as 

AF < B 2 ^spot _ B 2 1 / A phot AF mt \ 3 / 2 3/2 3/2 

A-^flare < — « — I j= j= rr 1 OC A t (X (A^ r ot) • (11) 

87T -L star T 87T L star r VI - (J spot/ 1 phot J ' 
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This indicates that the upper limit of the flare amplitude is proportional to 1.5 power of 
the starspots area if we assume B = constant. We discuss in the following whether this 
assumption is right. 

The discussion described around equations (ITlWffTTT) does not take into account the 
effect of the inclination angle i and spot latitude. In stars with lower inclination angle 
and higher spot latitude, the amplitude of brightness variation is smaller. If we consider 
inclination effects and assume the spot is distributed around the equator, equation ffTTj) can 
be replaced by equation f[T2"j) . 

B 2 (A spot / sin*) 3 / 2 B 2 1 / A phot AF rot 1 n^ 



AF flare < ^~»^< —"> ~r : t jii " L — IUL _ . (12) 

87T L star r 87rL star r VI- (T spot /T phot ) 4 sini/ 

If we adopt L star = 10 33 erg s" 1 , r = 10 4 sec, T spot = 4, 000K and T phot = 6,000K for 

superflares on solar-type stars, then the flare amplitude can be roughly estimated by 

B \ 2 / A 1 \ 3 / 2 



af «- k10 (iooog)Kw) ' (13) 



Figure [9] is the scatter plot of the superflare amplitude as a function of the ampli- 
tude of the brightn ess variation. The data of superflares are taken from quarter 1-6 data 



and are reported in 



Shibayama et al. 



(120131 ). The lines in this figure represent the analytic 
relation between the flare amplitude and the amplitude of the brightness variation obtained 
from equation ( Ti~3l) for i=90 degree and 2=2 degree (nearly pole-on) case each. Tak ing into 



accou nt that magnetic energy density of solar-type stars is about 1, 000 ~ 4, 000G (ISolanki 



20031 ). we show two lines (_B=1,000G & 3,000G) each. These lines give an upper limit 
for the flare amplitude based on the above discussion. Many of the flares have smaller 
amplitudes than that expected by the equation ffl3l) . Moreover all flares have smaller flare 
amplitudes than that expected in the nearly pole-on case (dashed-line; i—2 degree). This 
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supports that the energy released by these flares can be explained by the magnetic energy 
stored near the starspots. 



For the solar flares, there is a similar relation as this (ISammis et all 120001 1 . Figure 



[TUl sho ws the empirical re 



flares (ISammis et al. 



ation between the spot group area and X-ray intensity of solar 
20001 ). and the relation between the spot group area (estimated from 
the brightness variation amplitude) and the superflare energy. It should be noted that 
while there is a large range in observed flare energies for a given sunspot area, we here in 
particular refer to the upper envelope of the observations shown in Figure [TUl where the 
largest flare energies are associated with larger spots. The bolometric luminosity and the 
total bolometric energy of superflares were estimated from th e stellar radius, the effective 



temperature in the Kepler Input Catalog (IBrown et al. 



20111 ). observed amplitude and 



duration of flares by assuming that the spectrum of white light flares can be descri 



jed b y 



the blackbody radiation with the effective temperature of ~10,000K (IKretzschmar 
Using concrete numerical values, equation (jHJ) is transformed as 

B 2 



201lh . 



K 



in, « 



-.4 



3/2 



spot 

W A [erg\ 



in 
-,33 r 



10 33 [er#] 



By, A spot y /2 

10 3 G/ V3 x 10 19 cm 2 y 
B \Y ^spot/(27rg| )\3/2 

10 3 G7 V 0.001 



> E, 



flare 



(14) 



where R Q is the radius of the Sun. The relation expressed by equation f JT4|) can be seen 
in Figure [10J All of the solar flare data, and over a half of the superflare data are located 
below the same line of I?=1,000G, and i=90deg. 



Some data points of superflares are located above the line of equation (THj) . This is 
explained by the two reasons discussed in the following. First, the inclination and latitude 
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of the spots govern their projected area, and the size of the spot cannot be uniquely 
determined because of these effects and the fixed contrast with the photosphere. Lower 
inclinations cause a smaller projected spot area and/or the spot to be visible through most 
of the rotation period, and therefore lower brightness variations result in these cases. 

Second, magn etic flux densit y (B) is not necessarily constant. It is a function of 



the spot size. 



Solankil ( 120031 ) shows that the umbral continuum intensity of sunspots 



becomes weak as sunspots become large and that the intensity also becomes weak as 
the magnetic flux density gets strong. This fact means that as the spot becomes larger, 
the magnetic flux density becomes strong and the magnetic energy stored near the spot 
increases. As a result, the released flare energy tends to be bigger than that in the case of 
B=constant, and the straight lines in Figure [TD] become steeper. 

In this paper, we showed that many of superflare-generating stars probably have 
large starspots, and the energy of superflares can be explained by the magnetic energy 
stored around such starspots. We also confirmed in Section H] the stellar activity depends 
on the rotation period. Stars with relatively slower rotation rates can still produce flares 
that are as energetic as those of more rapidly rotating stars, although the average flare 
frequency is lower for more slowly rotating stars. There are, however, many things to be left 
for consideration in the next step. For example, it is extremely important to know lifetime 
of starspot regions or how activity level cha nges because it is rela ted to knowing how the 



starspot evolves by the dynamo mechanism. 



Shibata et al. 



(120131 ) indicated that the mag- 



netic energy to generate superflares can be produced by the effects of differential rotation. 
We need to investigate in detail the effects of differential rotation on the brightness variation. 



In addition, if superflares occur on the starspot area of a single star and the bright- 
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ness variation is caused by the rotation, the timings when superflares occur are likely to 
have a phase dependence, since superflares which we can observe need to occur when the 
starspots are on the visible side of the star. This phase dependence, however, is expected 
not to be necessarily true. For some pairs of spot distribution on the stellar surface and 
stellar inclination angle, some part of the area covered by the starspots is always on the 
visible side of the star. In these cases, it is probable that superflare occurs in any phase 
of the differential rotation . We need to study these effects in detail in the next step (cf. 



Roettenbacher et al 



2013) 



6. Summary 

In this paper, we investigated the brightness variations of superflare-generating stars in 
Kepler quarter 0-6 data by considering the signature of stellar rotation and the starspots. 
First, we performed simple spot model analyses for typical superflare-generating stars. 
Many of the brightness variations with the period of one day to a few tens of days and 
the amplitude of 0.1 ~ 10% can be explained by the rotation of the star with fairly large 
starspots, taking into account the effects of the inclination angle and the spot latitude. 
Next, using the result of the period analysis, we investigated the relations between the 
energy and frequency of superflares and the rotation period by assuming that the period 
of the brightness variation corresponds to the rotation period. Stars with relatively slower 
rotation rates can still produce flares that are as energetic as those of more rapidly rotating 
stars, although the average flare frequency is lower for more slowly rotating stars. Last, 
we discussed the relations between the starspot coverage and the energy of flares. If we 
assume that the brightness variations are due to the rotation, the energy of superflares can 
be explained by the magnetic energy stored near the large starspots. 
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Table 1: Length of the observation period during each quarter and the number of solar- type 
stars. 



Quarter 


N a 


r[days] b 


Release 





9511 


11 


14 d 


1 


75598 


33 


14 d 


2 


82811 


90 


14 d 


3 


82586 


90 


14 d 


4 


89188 


90 


14 d 


5 


86248 


95 


16 e 


6 


82052 


90 


16 e 



a Number of solar-type stars. 

b Length of the observation period during each quarter. 

c Release number of the data we use in this paper. 



d Kepler Data Release 14 Notes (JChristiansen fc Barclay 



Kepler Data Release 16 Notes (JChristiansen fc Barclay 



2012a ) 



2012bl ) 
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Table 2: Stellar parameters for three solar-type superflare-generating stars of which we show 
the results of model analysis in Figure [3j H] & [5j 



Kepler ID 


rp a 

J-eS 


log# a 


R/R & & 


Period b 


Amplitude 


6034120 


5407 


4.7 


0.77 


5.6 


4.45xl(T 3 


6691930 


5348 


4.5 


0.95 


13.1 


1.29xl0~ 2 


10528043 


5143 


4.5 


0.88 


12.9 


6.52xl0~ 3 



201 lh . 



a These data are taken from the Kepler Input Catalog (jBrown et al. 

fc The period of the largest amplitude component of the brightness variations in unit of day. 

c Thc normalized amplitude of the dominant period. 



Table 3: Best set of model parameters for Figure El @] & EJ 



Kepler ID 


Inclination 


spot name 


•f^spot / -f'phot 


Initial Latitude 


Initial Longitude Difference 


6034120 


45° 


spotl 


0.13 


41°N 


- 


6691930 


32° 


spotl 
spot 2 


0.23 
0.21 


30°N 
60°N 


70° 


10528043 


60° 


spotl 
spot 2 


0.15 
0.10 


0°N 

37°N 


120° 
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0.05- 




-0.05 



50 100 

Modified Julian Days - 54954.5 



Fig. 1. — Light curve of typical superflares on solar-type stars (KIC6034120 & KIC6691930, 
Quarter 2; KIC10528093, Quarter 0-2). The vertical axis is the brightness variations relative 
to the average brightn ess. The typica l photo metric error is about 0.02 %. Arrows indicate 



superflares reported in 



Maehara et al 



(2012). 
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(a) KIC60341 20 spectrum 
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(c) KIC1 0528093 spectrum 




frequency [day ] 



Fig. 2. — Power spectra of the light curves shown in Figure [TJ The periods of these three 
stars are shown in Table [2j 
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Fig. 3. — (a): Model light curve for KIC6034120 (Figured] (a)). The model parameters are 
given in Table |3j (b) & (c): Model pictures of the visible area of the photosphere with a 
starspot. (d): Observed light curve (solid line; the same as in Figured] (a)) and the model 
one (dashed line) for KIC6034120. 
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Fig. 4. — (a) : Model light curve for KIC6691930. The model parameters are given in Table 
[3j (b) & (c): Model pictures of the visible area of the photosphere with two starspots. (d): 
Observed light (solid line; the same as in Figured] (b)) and the model one (dashed line) for 
KIC6691930. 
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Fig. 5. — (a) Model light curve for KIC10528093. The model parameters are given in Table 
[3j (b), (c), (d) & (e): Model pictures of the visible area of the photosphere with starspots. 
(f): Observed light curve (solid line; the same as in Figure [TJ(c)) and the model one (dashed 
line) for KIC10528093. 
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Fig. 6. — Power spectra of the period analysis for the observed light curve and for the model 
one of KIC10528093 shown in Figure 0(f) . 
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Period (days) 



Fig. 7. — (a) Scatter plot of the flare-energy vs the brightness variation period. The period of 
the brightness variation in this figure was estimated by using the Kepler data of Quarter 1-6. 
We do not use Quarter data here since the data period is shorter (~11 days) than other 
Quarters and sometimes long term variations (~30 days) are removed. An apparent negative 
correlation between the variation period and the lower limit of the flare energy results from 
the detection limit of our flare-search method, (b) Distribution of the occurrence of flares in 
each period bin as a function of the brightness variation period. The vertical axis indicates 
the number of flares with energy > 5 x 10 34 erg per star per year. The error bars represent 
the la uncertainty estimated from the square root of the event number of flares in each bin. 
The frequency distribution of superflares saturates for periods shorter than a few days. A 
similar saturatio n is observed for the relations between the coronal X-ray activity and the 



rotation period (JRandich 



2000). 
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Fig. 8. — (a): Normalized irradiance of the Sun and area not coverd by sunspots. The 
data period is 2004.7.1 - 2005.6.30. Solid lines are normalized visible (4,500-8,000A) solar 
irradiance from the solar spectral irradiance data of the SORCE Satellite. We use the 
daily sunspots area data prepared by the U.S. Dept. of Commerce, NOAA, Space Weather 
Prediction Center, (b): Normalized irradiance of the Sun and the luminosity estimated from 
the spot coverages assuming that the temperature of sunspots is 4,000K and that of the solar 
photosphere is 6,000K. The data period is the same as (a). 
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Fig. 9. — Scatter plot of superflare amplitude as a function of the amplitude of the bright- 
ne ss variation. The data of superflares are taken from quarter 1-6 data and are reported 

in 



Shibayama et al. 



( 120131 ). We here defined the amplitude as the normalized brightness 
range, in which the lower 99 percent of the distribution of brightness difference from the 
average, except for the flares, are included. Thick and thin solid lines correspond to the 
analytic relation between the stellar brightness variation amplitude (corresponding to the 
spot area) and flare amplitude (correspond to the flare energy) obtained from equation (1131) 
for S=3,000G and 1,000G. The thick and thin dashed lines correspond to the same relation 
in case of nearly pole-on (i=2.0 deg) for £>=3,000G and 1,000G. These lines are considered 
to give an upper limit for the flare amplitude. 
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Fig. 10. — Flare energ y vs spot coverage of superfla res on solar-type stars (fi 



Shibayama et al 



20131 ) and solar flares (filled circles; 



Sammis et al 



2000 



Ishii 



ed squares; 



20121 ). The 



data o 



solar flares are take n from ftp://ftp.ngdc.noa a.gov 



1997 (jSammis et al. 



2000) and those in 1996-2006 fllshii 



/STP , and consists of data in 1989- 



20121 ). Thick and thin solid lines 



corresponds to the analytic relation between the flare energy obtained from equation ([14"]) for 
_B=3,000G and 1,000G. The thick and thin dashed lines correspond to the same relation in 
case of nearly pole-on (i=2.0 deg) for 5=3,000G and 1,000G. These lines are considered to 
give an upper limit for the flare energy (i.e., possible maximum magnetic energy which can 
be stored near sunspots). Note that the superflare on solar-type stars is observed only with 
visible light and the total energy is estimated from such visible light data. Hence the X-ray 
intensity in the right hand vertical axis is not based on actual observations. The energy of 
solar flares is based on the assumption that the energy of XlO-class flare is 10 32 erg, X-class 
10 31 erg, M-class 10 30 erg, and C-class 10 29 erg, c onsidering previous observational estimate 



of energies of typical solar flares (e.g., 



Benz 



2010l ). The values on the horizontal axis at the 



top show the total magnetic flux of spot corresponding to the area on the horizontal axis at 
the bottom when 5=1,0006. 



